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We present FaRe, a package for Mathematica that implements the decomposition of a generic 
tensor Feynman integral, with arbitrary loop number, into scalar integrals in higher dimension. In 
order for FaRe to work, the package FeynCalc is needed, so that the tensor strncture of the different 
contribntions is preserved and the obtained scalar integrals are gronped accordingly. FaRe can prove 
particnlarly useful when it is preferable to handle Feynman integrals with free Lorentz indices and 
tensor reduction of high-order integrals is needed. This can then be achieved with several powerful 
existing tools. 
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I. INTRODUCTION 

The computation of Feynman integrals is a very active research field and many tools are able to perform reduction 
and computation (see for instance m)- Very often, these routines require numerators with loop momenta contracted 
with outer momenta. While this is not an obstacle to any integral computation, there may be cases in which explicit 
Lorentz indices are preferable or, even more in general, the tensor decomposition of a high-order tensor Feynman 
integral is required. Different ways to reduce a tensor integral to easier, known scalar ones are available, such as the 
well-known Passarino-Veltman reduction [^, implemented, for example, in LoopTools m- Nevertheless, the higher 
the order and the larger the number of outer legs, the harder the decomposition becomes. A reduction technique 
proposed in Refs. HMS] appears very powerful and sufficiently machine-friendly to be implemented algorithmically. 
The package we propose, FaRe, deals with Feynman integrals, with arbitrary number of loops. The propagators 
can appear to any power and the tensor integral to decompose can be any rank. The output generated will be a 
combination of scalar integrals, in higher dimension, with different powers of the propagators. Most importantly, the 
tensor structure is preserved, so that it is possible to identify all the different Lorentz contributions in the final result 
and use it for further calculations in Mathematica. 

FaRe requires the package FeynCalc |16j in order to perform the decomposition and supply the result. By means 
of a built-in function already implemented, it is possible to get an output suitable for further reduction to master 
integrals with FIRE [7]. It is the usage of FaRe in combination with other tools, such as FIRE, that can prove to be 
rather useful in many different ways, even though the results provided by FaRe alone are always self-consistent. 

This paper is structured as follows: in Section|ll]we give a brief overview of the algorithm, in Section [llT] we describe 
the package FaRe and in Section |IV| we draw the conclusions. Finally, in the Appendix an explanatory example is 
considered. 


II. THE ALGORITHM 


FaRe implements the tensor reduction algorithm proposed by mHE], that makes use of raising and lowering 
operators. The procedure is briefly reviewed in what follows. 

As a starting point, we consider the iV-loop scalar Feynman integral in D-dimensional Minkowski space 


tD 


[fr d^ki 1 


( 1 ) 


where ki are the loop momenta and is a propagator whose generic exponent is I'i. This integral can be re-written 
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using the Schwinger parameterisation as 


tD 


poo " \ / " 


. ffr d^ki 

J *7r^/2 


where, in general, we have 
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1^1 h^l 


Z=1 


Here a; and cih are polynomial in Xi, while df and / are functions of both Xi and the outer momenta. This expression 
can then be diagonalised by performing a momenta transformation. It is convenient to transform the loop momenta 
in sequence. Calling C^{ki) and C{kf) the coefficients of the Z-th momentum and of its square respectively, we can 
replace k^ with the new momentum 


' “ ' 2C{kf)- 


( 4 ) 


It is clear that the first momentum to be transformed will have the trivial values C^{ki) = 2d'^ + 2J2j=i cij kj and 
C(kf) = oi, but the following replacements will have a more complicated structure. Once all the N replacements have 
been performed, the exponent will be 


N 


^^XiAi — '^^C{k‘f)Ki + — 


i=l 


l=l 


( 5 ) 


where Q/P regroups all the terms zero-order in all the loop momenta. The integration over the new loop momenta 
Kds can then be performed as standard gaussian integration, so that we end with 
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, x'''~^dxi 


= Vxl, 


i=l 
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n^^) * 


pD/2 


exp 
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( 6 ) 


where we defined 


Vx = 


rA(-ir' ..- 1 . 


( 7 ) 


We can identify a tensor integral according to the momenta that appear in the numerator, by giving the list M of 
their labels. For instance, a tensor integral like 




d^ki fc'' kl k^ 




( 8 ) 


can be specified by V = (1,1, 5), so that we can refer to the elements of the list as Afj (e.g. Afi = 1, etc.). Clearly, 
the order in Af is arbitrary. From eq. ([^ and ^ , a generic tensor integral of rank TZ can be obtained as 


«. ..(V-™) = /n^ 

1=1 j=l \i=l 


where we defined 


^ 

~ d{d^^ 




Performing on eq. § the diagonalisation and gaussian integration discussed above, we get 
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The action of the derivative operations on the exponential will result in a function which contains all 

the allowed tensor structures that can be formed with the external momenta Pe, and whose coefficients are functions 
of the Schwinger parameters Xi. We have 

= J ( 12 ) 

All the coefficients in W^^"'^'^{xi]Pe), polynomials in Xi, can be absorbed in the integration measure. Specifically, 
each Xi gives a factor —Vi and raises the relative exponent vi ^ Vi + 1 , while each factor 1 /P modifies the dimension 
of the integral D ^ D + 2 . Therefore, the tensor integral we started with is decomposed in a combination of scalar 
integrals in higher dimension, with different propagator exponents and regrouped according to the tensor structures 
given by the outer momenta. 

For detailed examples, we refer the reader to the literature (e.g. US]). 

FaRe automatically implements this procedure with arbitrary loop number N and tensor rank TZ. 


III. THE PACKAGE 


The package FaRe can be downloaded at https : //sourceforge .net/projects/f eyntoolfare/. It does not require 
any installation and it can be called in your Mathematica notebook simply by typing <<FaRe .m. It requires the package 
FeynCalc, which must be loaded in the notebook, lest FaRe loading automatically aborts. 

The functions available in FaRe and their arguments are listed below. 

1 
2 
3 


QP [ loops ] 

TRed [ D,num,den,loopMomenta ] 
FIREType[ expr,zeros,irr,fLetter ] 


The function QP returns the form ol Q/P in eq. (II), after the momentum transformations in eq. Q. The quantity 
returned does not contain the zero-order term / (see eq. (§), which can be easily added by hand. The output is in 


terms of coefficients d^* and 


where 


and 


= Cij for i ^ j. The only argument of the QP function is 


the number of loops. Here we give as an example the form of Q/P in the two-loop case 


QP[ 2 ] 


(acll d2^ -2 acl2 dl d2 +ac22 dl^)/(acll ac22 -acl2'') 


where we can see that P = 0102 — 0^2 and Q = aid/^ — 2 ci 2 di • ^2 + 02^1- Adding the / term, we get Q = 
aid| — 2 ci 2 dl • ^2 + a2d1 + /(ci2 — 0102), that matches the results found in the literature. 

The function TRed is the core of the package FaRe. It performs the tensor reduction of the integral provided. TRed 
has four arguments 


TRed[D,num,den,loopMomenta] 


• D is an integer representing the dimension D in which the integral must be computed. Notice that, using 
dimensional regularisation, the regulator must be left implicit and not specified in D. 

• num is a list containing the loop momenta that appear in the numerator of the integral. 

• den is a list containing the arguments Ai of the propagators A^, and their relative exponent in the form 
{{Ai,vi},. . . ,{An,Vn}}- Each Ai must be written without the square. 

• loopMomenta is a list containing the loop momenta. 

The output is a combination of scalar integrals gathered by the tensor structures formed with the external momenta. 
Integration over the Schwinger parameters is understood. The tensor structures that appear in the output are 
FeynCalc objects that can be manipulated in further computations. Indices and vector names are Mathematica 
strings. 

To illustrate how this function works, we can compute a double-box vector integral in three dimensions. Taking 
the momenta as in fig. we can compute the vector integral 
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(1,14.144.1) 1 1 y i7ri5/2 j7r^/2 ^2 1';,^ _ ^,2^2 ^2 

_ r d°ki d^k 2 _ fc^ __ 

J i7r^/2 iTT^/^ (fc^ _ k^f (fci - Pi2f (fcl - Plf kl {k2 + Pif {k2 + p^f ’ 


(13) 

(14) 


where pi2 = pi +P2 and P34 = p^ +pa. In the second step we have cancelled a propagator against k^ in the numerator. 
The corresponding FaRe input is 

1 
2 


Here, only the vector arguments of the propagators appear in the list, together with the relative exponent, in this 
case all equal to 1. The outer momenta are automatically converted into FeynCalc four-vectors by FaRe and treated 
accordingly. The output produced by FaRe is 


dbl=TRed[3, fkl}, {{kl-k2,1},{kl-pl2,1},{kl-pi,1>, 
{k2,1},{k2+p4,1},{k2+p34,1}}, {kl,k2}] 



where both the momentum name and its index are strings. 

The output produced by TRed can be further manipulated in Mathematica, and it can be used as input to FIRE, in 
order to reduce the obtained scalar integrals to master ones. To this aim, the function FIREType can be employed. 
This function has 4 arguments 

FIREType[expr,irr,zeros,fLetter] 


• expr is a linear combination of scalar integrals as obtained from TRed output. It is convenient to remove the 
tensor structures, even if it is not strictly necessary. 

• irr is the number of irreducible numerators treated by FIRE. Their exponents will be appended to the ViS lists 
and will appear as zeros. 

• zeros is the list of positions where to insert a zero in the t'^’s list. This can be used if FIRE asks for a list of 
exponents longer than that used by FaRe, for instance because of simplifications. 

• fLetter specifies the name used to represent the integrals in a suitable way for FIRE. 
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Referring to the previous example, we can apply FIREType to the coefficient of . FIRE will accept the seven 
denominators as in eq. (13) plus two irreducible numerators, so that it will handle lists with nine elements. Due to 
the simplification in eq. (14) the list produced by TRed have only six indices while, at the same time, TRed assures we 
do not have numerators anymore. The input and output of FIREType are then 


5 

6 

7 

8 
9 

10 


cpl = Coefficient [dbl , Pair [Lorentzindex ["/rl"] , Momentum [ " pi " ] ] 

~ (^{1.1,1,1,2,2} ■*" ^{1.1,1,2,1,2} ■*" ^{l.l,2,l,l,2}) 

FIREType[cpl,2,{3},G] 

-G [{1 , 1,0,1 , 1,2,2,0,0}] - G [{1 ,1 ,0,1,2,1 ,2,0,0}] 

-G [{1 , 1,0,2,1 , 1,2,0,0}] 


where two more indices are added and a zero is inserted at position 3, due 
arguments we have specified that the integrals will have function name G. 
once converted to expression. It must be noticed that FIREType does not 
be separately taken care of. 

An explanatory example on the usage of the functions available in FaRe 
with FeynCalc and FIRE, can be found in the Appendix. 


to the cancellation in eq. (14). In the last 
The output can be used as input to FIRE 
specify the integral dimension, that must 


and, more specifically, on its employment 


IV. CONCLUSIONS 


We have presented the new package for Mathematical FaRe, that automatises dimensional tensor reduction of 
an arbitrary-rank tensor Feynman integral with arbitrary number of loops. Its usage can prove to be useful for 
computation of high-order tensor Feynman diagrams, providing a rather easy way to obtain their tensor decomposition. 
The subsequent reduction of the resulting scalar integrals to master ones can be achieved by means of other powerful 
existing tools such as FIRE. The possibility of interfacing FaRe with FIRE is made easier through the function FIREType. 
However, improvements can be made in this respect. 

The synergetic employment of FaRe, FIRE and, of course, FeynCalc can thus prove to be quite powerful in handling 
high-order Feynman integrals. In the Appendix we sketch a possible example, borrowed from ABJM theory |I8j . 
that can be find in its detailed version in the example notebook attached to the FaRe distribution. Here we compute 
the a double-box integral in three dimensions. Though already well-known, the result obtained with FaRe does not 
represent just a benchmark of the code, but shows hidden identities which may prove to be interesting (see for instance 
eq. (AI5|). 

The package FaRe is fully open source, which represents one of its strong point. Open access to the source code 
allows for improvements from the whole scientific community. This way, the door stays always open for countless 
customisations that may help to tackle many different problems in high energy Physics. 
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Appendix A: A pedagogical example: double-box in three dimensions 

As an explanatory example of the functioning of FaRe and its interplay with other reduction tools such as FIRE, 
in this appendix we compute a double-box integral entering the planar two-loop ABJM amplitude [IMI]- The full 
computation is reported in the example Mathematica notebook that can be found together with the package FaRe. 










6 


Defining the momenta as in fig. we compute the integral 

d^ki d^k2 


DBp(s, t) = 


Afp 


It,D/1 j^C/2 ^ ^2 fc2 (k2 + PSA? (fe + Pif (fcl - k2f (fci - pif ’ 


(Al) 


where 


AAp = [str{fcipip4} + fcj tr{pip2P4}] [str{fc2PiP4} + fci tr{piP2P4}] ■ 


(A2) 


The integrals to be decomposed by FaRe are just three, since the product of the last two terms in both traces gives 
a scalar integral with no loop momenta in the numerator (thanks to the cancellation against the two corresponding 
denominators). We have, therefore 


tr {fci pipi} tr {k 2 PiPi} 


f\ = — f _ 

^ ^ J kf {ki - pi2f k‘^{k2 + P34f ik2 +pa^ (ki - k2f {ki - pif’ 




DB^^^(s,t) = 


■D/2 i-^D/2 ^2 ^ 

str{piP2P4} f d^ki d^k2 


tr{k2PiP4} 


t 


i 7 r ^/2 mD/2 ^2 p ^^2 _ ^^^2 _ ^^^2 


str{pip2P4} f d^ki d^k2 


tr{fcipip4} 


t 


i^D/2 ij,D/2 j^2 ^ ^^^^2 ^^^2 _ ^^^2 _ ^^^2 ’ 


(A3) 

(A4) 

(A5) 


where in eqs. (A4) and 


we have performed cancellations between loop momenta in the numerator and in the 


denominator. Recalling that the trace yields a Levi-Civita tensor, we have to compute 

‘'d^ki d^k2 k^k)^ 


hi)(Kk2) = 


1 ( 2 ) 


- /) 


d^ki d^k2 


k^2 




d^ki d^k2 


K 


iT,D/2 _ p^^^2 ^2 ^^^2 _ ^^^2 


\2 ’ 

(A6) 

pi) 

(A7) 

r 

\2 ’ 

(AS) 


of which only terms proportional to p 2 , Pa and the metric tensor survive. Solving the momentum conservation relation 
for p 3 , we can just keep terms proportional to g^'' and P 2 - The integrals in eqs. (A7) and (ASI are vector integrals 
already computed in Section [iTj They can be easily decomposed by FaRe as follows 


12 = 

TRed[3,{k2},{{kl-pl2,l},{k2, 

1>,{k2+p34,1},{k2+ 

p4, 1 },{kl- 

-k 2 ,l}. 



{kl-pl, 1 }},{kl,k 2 }]; 






f 12 

= 12 /. {Pair[Lorentzindex[ 

x_] , 

Momentum["p34 

" 1 ] : > 




FourVector["p3" 

. x] 

+ F ourVector[ 

"p4", x], 




Pair[Lorentzindex[ 

x_] , 

Momentum [ "p12 

" 1 ] : > 




FourVector["pi" 

, x] 

+ FourVector [ 

"p 2 ", x]} 

// Fulls 

implify ; 

FI2 

= fI2 /. Pair[Lorentzindex[ 

x_] , 

Momentum [ "p3 " 

]] :> 




-FourVector["pi 

" , x] 

-FourVector [ " 

p2", x]-F 

3 urVe ct or 

["p4" , x] ; 

a 2 

= Coefficlent[Expand[FI2], 







Pair[Lorentzindex ["ml"] ,M 

oment 

um["p 2 "]]] // 

FullSimplify 



—T® _ 7-5 4- T® 

{ 1 , 1 . 2 , 1 . 1 , 2 } -^{ 1 , 1 . 2 , 1 . 2 , 1 } ^-^{2 

, 1 . 1 . 1 - 2 , 

L} -^{2,1.2,1.1,1} 




A2 = 

FIREType[a2,2 , {1} , G] // ToEx 

press 

ion 





-G[{0,1,1,2,1,1,2,0,0}] 

-G[{0 

112 12 1 

0 , 0 }] 




+G[{0,2,1,1,1,2,1,0,0>] 

-G[{0 

2 12 111 

0 ,0}] 




I3=TRed[3,{kl},{{kl,l}{kl-pl2,1>,{k2+p34,1},{k2+p4,1},{kl-k2,1}, 
{kl-pl,1}},{kl,k2}]; 

fI3 = 13 /. {Pair[Lorentzindex[x_], Momentum["p34"]]:> 

FourVector ["p3", x] + FourVector["p4", x] , 

Pair[Lorentzindex[x_], Momentum["pl2"]] :> 

FourVector["pi", x] + FourVector["p2", x]} // FullSimplify; 
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FIS = fI3 /. Pair[Lorentzindex[x_] , Momentum ["p3"]] ;> 

-FourVector [ "pi", x]-FourVector["p2", x]-FourVector["p4", x] ; 

a3 = Coefficient [Expand[FIS] , 

Pair [LorentzIndex["ml"],Momentum["p2"]]] // FullSimplify 

-^{1,1,24,2,1} ~-^{1,2,1,1,2,1} ~-^{1,2,1,2,1,1} ~-^{1,2,2,1,1,1} 

A3=FIREType[a3,2,{3},G] // ToExpression 

G[{1,1,0,2,1,2,1,0,0}]-G[{1,2,0,1,1,2,1,0,0}] 

-G[{1,2,0,1,2,1,1,0,0>]-G[{1,2,0,2,1,1,1,0,0}] 


The computation of /(i), that is a rank-2 tensor integral, is less trivial. Thanks to what discussed above, we are 
only interested in the terms proportional to the metric tensor and p 2 - We have 


Il=TRed[3,{kl,k2>,{{kl,l},{kl-pl2,l},{k2,l},{k2+p34,1},{k2+p4,1}, 

{kl-k2,1},{kl-pl,1}},{kl,k2}]; 
fll = II /. {Pair[Lorentzindex [x_] , Momentum["p34"]]:> 

FourVector["p3" , x] + FourVector["p4", x] , 

Pair[Lorentzindex [x_] , Momentum["pl2"]] :> 

FourVector["pi", x] + FourVector["p2", x]} // FullSimplify; 

Fll = fll /. Pair[Lorentzindex [x_] , Momentum["p3"]] ;> 

-FourVector["pi" , x]-FourVector["p2" , x]-FourVector["p4" , x] ; 

al=Coefficient[Expand[FI1], MetricTensor["ml", "m2"]] // FullSimplify 

7-5 

-^{1,1,1,1,1,2,1} 

Al = FIREType[al, 0, 2,G] // ToExpression 

-((G[{l,l,l,l,l,2,l,0,0}])/2) 

(* dimension = 7 here *) 

bl = Coefficient [Expand[FI1] , Pair[LorentzIndex["ml"],Momentum["p2"]]Pair [ 
Lorentzindex["m2"],Momentum["p2"]]] // FullSimplify; 

Bl=FIREType[bl,2,{0},G]; 

The expression for Bi is rather lengthy and can be found in the example notebook. Nonetheless, it must be noticed 
that Bl is the only coefficient whose integrals are in dimension 7, while for all the other quantities we have D = 5. 
It is now possibile to reduce the scalar integrals to master ones, using FIRE. The details of the decomposition are 
reported in the notebook. We show here the analytical expressions of i?i, A 2 , and A 4 


Bl = 

^2 = ^3 = 


27r r 1 
\ e 
27r r 1 
1 2 e 


s-2t(l-Flog4)-P 


s - t (1 -F log 4) -F 


2 t log t 


tlogt 


+ 0 {e) 
FO(e 


Ott 

A 4 = -^+ 0 {e). 


(A9) 

(AlO) 

(All) 


Here A 4 is the scalar integral obtained from the product of the two last terms in the traces of eq. (A2). As already 
mentioned, it does not require any tensor reduction. 

It is important to highlight here a subtle feature of the dimensional reduction (DRED) employed so far. The space- 
time dimension of the theory imply that the external momenta are 3-dimensional vectors. However, the loop-momenta 
are chosen to be H-dimensional vectors, where H = 3 — 2e, in order to regularise the loop integrals with e. This 
dimensional mismatch shows up in the computation of the term proportional to Ai. Here we have 




(A12) 


where e denotes the Levi-Civita tensors. As already mentioned, only the terms proportional to and P 2 P 2 in 
survive, so that 


(S, t) = 'iepapf-upaPlPl P 4 P 4 (^1 V'"'" + Bl pt^p^) ■ 


(A13) 














Here, the metric tensor stemming out of DRED is implicitly £)-dimensional, while the Levi-Civita tensors and the 
outer momenta are three-dimensional. The contraction of the two e tensors with the Z?-dimensional metric yields 
a factor (1 — 2 e) [52], so that the complete integral is then given by 

DBp(s,t) = (1 — 2^) s'^tAi + s^{s + t) {—sBi + A2 + A3 — A4) . (A14) 

Normalisation factors are absorbed into the master integrals in Hi_ 2 , 3,4 and Bi. It is interesting to notice that the 
second term vanishes up to terms of order e, i.e. 


— s Bi -j- A 2 A 3 — A 4 — 0(e), 


(A15) 


so that, eventually we get 


DBp(s, t) = (1 - 2e) s^tAi -b 0(e) 

= - 4 + [2(1 - log2)-l-log(st)]--b 


27r 


27r^ 

2 -b —- -b 2 log^ 2 -b log s (2 log 2 - 1 - log t) 

o 


0 {e), 


(A16) 


which matches the result found in the literature 
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